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I.  INTRODUCTION 


For  guided  missiles  targeted  against  certain  classes  of  targets,  underground  detonation  of  the 
warhead  may  be  more  effective  than  air  hurst  in  achieving  a desired  level  of  damage.  Also, 
underground  detonation  can  reduce  the  possibility  of  unacceptable  nonmilitary  casualties  and 
collateral  damage.  The  target  is  assumed  to  be  finite  in  the  sense  that  the  size  of  the  taiget  is 
small  in  relation  to  the  standard  deviation  of  the  distribution  of  miss  distances  possible.  If  the 
guidance  syst'^m  is  such  that  no  information  on  the  vector  miss  distance  between  the  weapon  and 
the  finite  target  is  available  or  can  be  derived,  no  target-oriented  warhead-detonation  logic  can 
be  prescribed.  Thus,  for  this  case  the  warhead-detonation  command  must  be  derived  from 
environmental  information  obtained  near  the  target.  For  example,  the  fuzing  logic  for  this 
application  could  be  to  command  detonation  at  a fixed  time  interval  after  earth  entry.  With  use 
of  integrating  accelerometers,  the  fuzing  logic  could  command  detonation  after  a prescribed 
travel  distance  or  at  a prescribed  depth  below  the  surface.  None  of  these  techniques  define  such 
(possibly  more  desirable)  detonation  points  a^  the  point  of  closest  approach  (PCA)  between 
missile  and  target  (for  isotropic  warheads)  or  within  an  acceptable  isodamage  contour  (for 
nonisotropic  warheads).  

However,  if  missile  guidance  parameters  with  respect  to  the  target  are  available  for  the  pre- 
entry portion  of  the  trajectory,  along  with  appropriate  sensors  for  measuring  trajectory  after 
earth  entry,  the  problem  of  fuzing  at  an  app'opriate  point  defined  with  respect  to  the  target 
becomes  deterministic.  It  is  th-^  purpose  of  this  study  to  examine  the  fuzing  requirements 
analytically,  and  to  derive  the  necessary  relationships  and  logic  for  fuzing  an  earth  penetrator 
weapon.  The  weapon  is  assumed  to  be  deployed  against  an  underground  target  or  against  a target 
where  below-the-surface  detonation  could  result  in  increased  effectiveness  or  in  a reduction  in 
collateral  damage,  while  still  achieving  a desired  level  of  kill.  For  the  purpose  of  this  study  we 
assume  some  form  of  terminal  guidance  that  is  target  oriented.  That  is,  an  estimate  of  the  vector 
miss  distance  b,:tween  the  missile  and  the  designated  target  is  available  or  can  be  derived  from 
data  which  are  also  used  for  terminal  guidance.  It  is  also  assumed  that  this  information  can  be 
provided  for  the  fuze  processor  at  the  time  the  measurements  are  made  or  within  a few 
milliseconds. 

Certain  factors  affecting  the  desired  performance  of  the  system  are  assumed  either  to  be 
known  before  launch  and  to  be  stored  in  the  fuze  memory,  or  to  be  provided  to  the  fuze  by  a 
communications  link  before  earth  entry.  One  such  factor  Meeting  the  system  performance  may 
be  a minimum  desired  burst  depth  below  the  surface,  to  satisfy  some  constraints  on  collateral 
damage,  that  is,  damage  to  possibly  nonmilitary  facilities  near  the  prescribed  target.  A second 
factor  assumed  to  be  known  to  the  fuze  is  the  location  of  the  target  with  respect  to  the  guidance 
aim  point.  For  example,  in  the  case  of  a target  completely  below  the  surface,  the  guidance  miss 
vector,  for  the  type  of  guidance  system  assumed  here,  would  be  measured  to  some  point  ot  the 
surface  of  the  earth.  The  location  of  the  desired  target  may  be  at  the  same  earth  coordinates 
(latitude  and  longitude)  as  the  guidance  aim  point,  and  at  a known  depth,  or  it  may  be  offset  by  a 
known  distance  and  direction  from  the  guidance  aim  point.  This  information  would  have  to  be 
made  available  to  the  fuze  in  the  form  necessary  to  permit  use  of  the  algorithms  based  on  the 
intercept  geometry  logic  described  below.  Another  factor  assumed  to  be  known  to  the  fuze 
beforehand  is  the  information  from  which  the  desired  burst-point  requirements  may  be 
determined  and  used  in  the  fuze  logic.  The  fuze  burst-point  logic  may  also  require  data  from  the 
fuze  sensor  to  function  as  desired. 

This  report  will  consider  several  cases,  depending  on  the  form’  of  the  underground 
trajectory,  the  type  of  target,  and  the  form  of  the  available  initial  data.  First  the  initial 


3 


parameters  as  derived  frum  guidance  data  are  discussed,  including  the  assumptions  as  to  the 
form  in  which  they  are  available  and  what  further  parameters  can  be  derived  from  them.  It  is 
assumed  in  general  that  the  missile  remains  underground  after  impact,  and  that  there  are  two 
k.  own  depths:  a minimum  desired  burst  depth,  m,  and  another  depth,  t/,  to  the  location  of  the 
defined  target  or  to  the  coordinates  of  the  center  of  some  isoviilnerable  volume.  The  time 
parameter  is  defined  to  be  zero  on  impact,  and  increasing  thereafter.  For  the  purpose  of  this 
analysis,  the  target  is  assumed  to  be  at  or  below  the  aim  point;  any  other  case  would  require 
adding  a coordinate  translation  to  the  derived  relationships. 

For  all  cases  considered,  a fuzing  logic  is  defined,  based  on  the  initial  conditions, 
underground  trajectory,  target  position,  and  collateral  damage  criteria.  Two  different  criteria  for 
target  destruction  are  assumed.  The  first  is  that  detonation  occurs  at  the  PCA  between  missile 
and  target.  The  second  is  that  detonation  occurs  within  some  lethal  contour,  defined  in  terms  of 
warhead  directivity,  and,  if  applicable  and  available,  soil  characteristics.  The  effect  on  the  logic 
of  imperfect  data,  imperfect  measurements,  and  incomplete  information  for  a deterministic 
solution  is  also  studied. 

In  the  analysis  that  follows,  several  simplified  targets  and  underground  trajectories  are 
considered.  The  intercept  geometry  is  first  analyzed  for  a linear  underground  path,  and  then 
generalized  to  a path  created  by  a generalized  round  acceleration.  The  target,  as  a first  step,  is 
considered  to  be  either  a point  or  a line.  Each  target  case  is  then  expanded.  The  point  target  is 
generalized  to  an  isoviilnerable  spherical  contour;  the  line  target  is  generalized  to  an  isovulncra- 
ble  cylindrical  contour. 


2,  INITIAL  CONDITIONS 

In  general,  the  coordinate  system  is  set  up  as  follows.  The  origin  is  on  the  designated  target 
for  the  point  target  case.  The  necessary  modifications  of  this  choice  for  linear  targets  and 
isovulnerable  spherical  and  cylindrical  contours  are  noted  where  relevant.  The  x„  axis  is  the 
surface  line  from  the  aim  point  to  the  impact  point.  Fpr  this  analysis,  the  surface  is  considered 
flat  in  a sufficiently  large  neighborhood  of  the  aim  point.  The  z axis  is  the  surface  normal  that 
passes  through  the  aim  point  and  the  target.  Next  the  y,  axis  is  chosen  so  that  x^-y'o-z  forms  a 
mutually  perpendicular  right-handed  triple.  In  this  fashion  the  v>o  plane  defines  the  surface 
plane  in  a neighborhood  of  the  aim  point.  Finally,  the  ah,  and  .Vo  axes,  when  translated  downward 
along  the  z axis  to  the  target  origin,  become  the  x and  y axes.  This  construction  is  diagrammed  in 
figure  I . 

With  respect  to  this  system,  the  necessary  parameters  for  this  analysis  can  be  defined  more 
precisely.  The  aim  point  is  the  point  (0,0,t/)  and  the  point  of  impact  is  {Rfi.d)  where  R is  the  miss 
distance  and  d is  the  depth  to  the  target.  The  vector  miss  distance  has  length  R and  is  directed 
from  the  aim  point  to  the  impact  point.  The  further  parameters  which  will  be  needed  for  this 
analysis  are  the  impact  angle  with  respect  to  the  surface  a,  the  map  coordinate  system  angles  cr 
and  y,  impact  veliKity  V,  deceleration  A,  and  the  surface  angles  6'  and  0. 

2.1  Parameters  0'  and  0 

The  value  of  O'  is  required  to  determine  whether  the  impact  is  early  or  late.  The  value  of  6 is 
needed  in  the  computations  of  the  optimal  time  for  detonation.  These  two  parameters  are  derived 
from  the  guidance  data  as  follows. 


/ 


/ 


The  missile  is  assumed  to  have  sensing  mechanisms  that  can  determine  the  angle,  y,  of  the 
miss  distance  vectors,  ,7,  with  respect  to  the  map  coordinate  system  of  the  aim  point.  Angle  y is 
defined  as  the  angle  measured  from  north  to  R counterclockwise,  and  is  assumed  to  be  always 
positive.  Further,  the  sensor  can  defirmine  the  angle,  <t,  of  the  surface  component,  5,  of  the 
trajectory  with  respect  to  the  map  coordinate  system.  Similarly,  «r  is  the  angle  from  north  to  S 
counterclockwise:  again,  o is  always  positive.  Angle  0'  is  defined  to  be  the  smaller  angle 
between  R and  S. 

Figure  2 shows  one  possible  relationship  between  <t,  y,  and  9'  (see  app  A for  diagrams 
illustrating  all  such  possible  relations).  Upon  examination  of  these  diagrams,  the  parameter  9'  is 
calculated  as  follows. 


!f0  s |<7' - y|  s ;r,  then  9'  = |o--y|,  or 

if  rr  < |cr  - y|  s 27r.  then  9"  = 2w  - |o- - y|. 

In  all  cases  0 s d'  k n. 

Next,  the  angle  9 is  defined  from  9*  by 

9 ■■  9'  if  0 ss  9*  s rr/2,  or 
9 * T — 9'  if  Tr/2  s 9' s JT. 


N 


Figure  2.  Initial  angle  parameters. 


(3) 

(4) 
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2.2  Karly  and  Late  Impacts 

To  define  criteria  for  early  and  late  impacts,  the  surface  projection  of  the  pre-entry 
trajectory  is  first  assumed  to  be  linear.  Let  L be  the  line  in  »hc  surfae.’  plane  through  the  aim 
pi>int  and  normal  to  the  surface  trajectory.  The  impact  is  defined  to  be  early  if  the  surface 
trajectory  dt>cs  not  cross  L before  impact.  For  example,  the  impact  point  of  path  1 in  figure  3 is 
early.  The  impact  is  defined  to  be  late  if  the  surface  trajectory  has  crossed  L.  The  impact  point  of 
path  2 in  figure  3,  for  example,  is  late.  If  the  impact  point  is  on  L,  the  missile  impact  is  defined  as 
late  also. 

If  the  surface  projection  of  the  pre-entry  trajectory  is  not  rectilinear,  the  tangent  line  to  the 
curve  at  the  point  of  impact  replaces  the  linear  surface  trajectory  above.  Now.  L is  the  normal 
through  the  aim  point  to  the  tangent  at  the  impact  point.  Early  and  late  impacts  are  then  defined 
as  for  a linear  surface  trajectory.  This  situation  is  illustrated  in  figure  4 for  the  case  when  the 
impact  is  late. 


Figure  3.  Impacjt  in  surface  plane:  (a)  early  versus  (b)  late. 


SURFACE  TRAJECTORY  COMPONENT 


Figure  4.  Curved  trajectory  impact  type. 


Impact” ’oint 
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2.3  Criteria  for  Impact  Form 


(5) 

(M 


If  It'  s n/2,  the  impact  is  late; 
if  O'  > n/2,  the  impact  is  early. 

The  details  of  the  derivation  of  equations  (5)  and  (6)  are  given  in  appendix  A.  Hiieny.  a set  of 
diagrams  similar  to  figure  2 are  analysed  for  impact  type  and  the  corresporuling  form  of  0' . 

2.4  Derivation  of  Mks  Distance 

In  the  calculation  of  the  miss  distance.  W.  it  is  assumed  that  the  impact  angle,  it.  between 
the  pre-entry  trajectory  and  the  surface,  as  well  as  the  angle  ; 0,  y.  and  care  known.  Then  either 
/t.  the  current  height,  or  «.  the  distance  along  the  ir:ijecti>ry  to  the  impact  point  at  any  particular 
time  before  entry,  is  required  to  determine  R.  Also  required  is  fi.  the  angle  between  the 
trajectory  and  the  line  to  the  aim  point  (see  fig.  5).  I'hese  parameters  arc  determined  at  .K'me  pre- 
impact time.  But.  since  determination  of  R also  requires  information  as  to  whether  in  patl  is 
early  or  late,  the  calculation  of  /?  must  be  done  after  impact.  1 he  two  possibilities  are  depi  :tetl  in 
figures  5 and  6.  The  segment/ is  constructed  to  be  normal  to  the  .v  axis.  From  these  figures,  the 
parameters  « or  li  satisfy 


<•'  " /i/sin  a. 

In  both  cases,  the  lengths  h.f.  and  g are  given  by 

(7) 

h - a cos  n 

(K) 

II 

(9) 

1,’  = h cos  0. 

(10) 

TRAJECTORY 


/ 


Figure  6.  Early  impact  geometry.  8 - n - 8'. 


2.4.1  Derivation  fur  Early  Impact 
Given 

fl«  - + /»>. 

* /,»  + e*, 

e*  + (R  + gf,  and 
R*  - a*  + c*  - 2ac  cos  j3. 

From  equations  (12)  and  (14). 

c*  ” h*  + f*  = + P + + 2Rg  + /?*. 

Using  equations  (II)  and  ( 13)  to  simplify. 

c*  « fl*  + 2Rg  + /?* 

Substituting  equation  (17)  into  equation  (13)  gives 

R*  a*  + (a*  + 2Rg  + R*)  - 2ac  cos  /3,  and 
ac  cos  p = («*  + Rg). 

By  equation  (17). 

a*  + Rg 

(a*  + 2Rg  + /?»)•'*  = c-  = ^ . 

acosp 

After  both  sides  are  squared,  this  simplifies  to  a quadratic  in  R, 

(a*cos*/3  - g*)R*  - (2ga*sin*p)R  - a*sm*p  = 0. 


(11) 

(12) 

(13) 

(14) 

(15) 

(16) 

(17) 

(18) 

(19)  j 

(20)  j 

i 

(21) 

I 

i 

I 
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Solving  this  quadratic  for  R and  simplifying. 


R = 


Ha 


* sin*  /3  sin  cos  fi  n/«*  - V* 


(rt*cos*/3  - V*)  (a*cos*/3  - A'*) 


2.4.2  Derivation  for  Late  Impact 

Given  equations  (II)  through  (13)  as  in  the  early  case. 

R = g + k 
= P + k^,  and 
A?*  = o*  + (•*  — 2ac  cos 
Then  from  equations  (i2),  (24),  and  (23). 

(.2  = /,2  + ^2  = 1,2  + J1  + 1^2  f,2  + p + 

Simplifying  bv  equations  (11)  and  (13). 

P = a*  + R^  - 2Rg. 

Eliminating  i * in  equation  (25), 

/?*  = <!*  + (a*  + /?*  - 2Rg)  - 2ac  cos  /3,  and 
ac  cos  p = P - Rg. 

Squaring  and  using  equation  (27), 

c*  = «*  - 2/?g  + A?*  = — . 

•rcos^P 

Rewriting  as  a quadratic  in  R, 

{a*cos*/3  - g*)A?*  + 2</*gsin*/3  A?  - <i^sin*/3  - 0. 
Solving  for  A?  as  before  and  simplifying, 

r gPsin^p  ^ a*  sin  ^ cos  ^ rp — 

L(a*cos*/3  - ^)  (a*cos*/3  - )?*)  J 

Combining  the  ee  ly  and  late  cases. 

a*  sin  )3  (g  sin  /3  ± cos  /3  - g^) 


R^-± 


(a*cos*/3  - p) 


where  the  (+)  value  of  the  first  (±)  sign  corresponds  to  an  early  impact  and  the  (- 
late  impact. 


From  figures  5 and  6, 


So  in  'erms  of  a,  0,  and  p. 


g = b cos  0, 
b = a cos  a. 

A?=  ±aK[L  ^ cosp  VmI 


(22) 

(23) 

(24) 

(25) 

(26) 

(27) 

(23) 

(29) 

(30) 

(31) 

(32) 

(33) 

) value  to  a 

(34) 


II 


where  a is  know'n  or « - h/sin  a if  /i  is  known, 

h'  = sin  /3/(cos*  /3  - cos*a  cos*  B), 

L = cos  n cos  (f  sin  fi,  and 
M - \ - cos*  et  cos*  B. 

The  initial  data  from  the  guidance  sensors  and  the  necessary  parameters  subscqui.ntly 
calculated  are  summari/ed  as  follows. 

I.  Necessary  logic  inputs 
a impact  angle 

tr.y  map  system  angles 

V impact  veliKity 

A deceleration  a.'tcr  entry 


miss  distance; 


(a)  given 

(b)  calculated  from  equation  (34) 


2.  Further  parameters  tit  be  calculated 

B",  the  smaller  positive  angle  between  R and  S 

(a)  S' = |<r  - yl.  if  |<r-ylsjr 

(b)  S'  -2n  - \<T-  y|,  if  2ir  ?■  |o-  - y|  > tr 

(c)  0 £ S'  s ff 

B,  the  acute  angle  defined  by 
ia)  B=  S'  if  S'  s n/2 
(b)  » = jr  - S'  if  irl2^B'  Sir 

3.  POINT  TARGET 

Using  the  pre -entry  data,  the  information  necessary  to  define  the  fuzing  requirements  can  be^ 
derived.  A minimum  burst  depth,  m,  and  the  values  for  a.  6,  y,  a-,  R,  J,  and  V are  assumed  to  be 
known. 

3.1  Case  1:  Linear  Underground  Path  to  Target 

For  a linear  underground  path  to  target,  a linear  path  both  above  and  below  ground  and  a 
constant  deceleration  are  assumed  (see  fig.  7).  The  parametric  equations  of  the  position  of  the 
missile  as  a function  of  time  are 

Jt(t)  = .*i,  + V,/  + y f»,  (35) 

y(0  = y.  + + y and  (36) 

z(t)  ^ + V,/  + y r*.  (37) 


I 


where  Aj..  A^,  are  the  components  of  the  deceleration  A,  where  t is  assumed  to  be  zero  at 


time  of  impact,  where  x„  = R,  Vo  = 0,  and  = d,  and  where 

V,  = V cos  a cos  6,  (38a) 

V„  = V cos  a sin  and  (38b) 

V,  = V'  sin  a,  (38c) 

and  similarly, 

A^  = A cos  « cos  8,  (39a) 

Ay  - A cos  a sin  0,  and  (39b) 

A,  = A sin  a (where  A <Q).  (39c) 


2 


Figure  7.  Linear  underground  trajectory  to  isovulnerable  spherical  volume  about  point  target. 

Next  the  specific  parametrizations  are  defined  considering  the  direction  of  change  in  each  of 
the  components  of  the  missile  position.  Since  z(r)  is  always  a decreasing  function, 

z(t)  * </  - V sin  o / + sin  (40) 

The  fUiiction  y{i)  is  of  the  form: 

y(t)  = ±V  cos  a sin  Ot  ±—  cos  a sin  9t*,  (41) 
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w here  Ihe  iippeT  sign  pertains  to  r'  inereasing,  and  the  low'cr  to  v decreasing.  The  function  .v(r)  is 
of  the  form 

.v(/)  — R + — coS(»  cos  l>l~  T V'coso  cos  lit.  (42) 


w here  the  upper  sign  pertain  < to  early  impaet,  and  the  lower  to  late  impact.  As  the  missile  moves 
tilong  its  underground  trajectory,  the  data  processor  must  determine  when  the  missile  will  be 
closest  to  the  target.  Further,  the  time  to  reach  the  minimum  burst  depth  must  be  determined. 
For  both  quantities  some  form  of  distance  function  must  Ik  used.  The  distance  from  the  target 
e.rigin  to  a given  point  on  the  subsurface  trajectory  is  given  by  /hr)  = V v"(/)  + y"(r)  + ;^(/).  To 
minimize  i>{t)  it  is  siiftlcient  to  minimize  />"(/),  by  setting  Ihe  time  derivative  equal  to  zero.  The 
time  at  which  the  missile  is  at  Ihe  FC'A  is  defined  as  t.  Computing  y~{t).  and  z'Ht)  from 

equations  (4(1)  through  (42)  gives  the  expression  for  /r(/)  in  equation  (43)  below.  Consequently, 
the  expression  for  is  not  affected  by  the  possible  increase  or  decrease  in  either  the  .r  or  y 
component.  However,  the  early  ;ind  late  impact  ctises  do  inlluence  Ihe  expression.  The  upper 
signs  pertain  to  an  early  impact  situation;  the  lower,  to  ti  late  impact  situation.  This  sign 
convention  will  be  follov^ed  throughout  the  discussion  unless  otherwise  noted. 


^"‘t)  - t*  + AVi  ''  + (F^  f /^A  cos  II  cos  a - AJsino:)/^ 


(43) 


+ ( ■*■  2 V7?cos  llcos  a ~ 2 Vd sin  a)t  + R‘  i-  d'^. 


Jf)“ 

Then  = 0 implies 


0 = A't'  + 3A  W'  + 2(  V*  + RA  cos  d cos  or  ~ Ad  sin  o)/ 


(44) 


+ (-2V’</sin  or  + 2VR  cos  or  cos  d). 


Dividing  substituting  .v  - V/A  for  /.  and  solving  for  .v  in  equation  (44)  gives 

.V  = 0 

or 

.V  = ± — VF'^  + 2A(±K  cos  (f  cos  d f </sin  o'). 

A 

Then  t must  be  chosen  from  one  of  these  values: 

/,  = — V'/ A + — \/F  + 2A{±R  cos  II  cos  a + d sin  or). 

A 

1.,  = -VIA.  or 

= - V'/A  — - \/V^‘  + 2A(±/^cos  llcos  a + d s\n  o). 

A 


(4.5) 


(45) 

(47) 

(48) 


For  each  of  these  three  values,  the  function  /F(/)  has  an  extreme  point.  The  value  for  /, 
represents  the  time  to  retieh  the  first  local  minimum  of  fr{l).  Next,  t-i  corresponds  to  the  time 
.’•equired  for  the  missile  to  re;ich  its  maximum  penetration  and  stop.  The  value  of  l,,  represents 
Ihe  improbable  situation  of  the  missile  having  negative  velocity,  retracing  the  path  out  of  the 
earth,  and  passing  through  the  local  minimum  again.  To  determine  when  to  detonate  the  warhead 
so  that  all  necessary  conditions  are  satisfied,  the  time,  r,  to  reach  the  closest  achievable  point  to 
the  target  must  first  be  calculated.  Then  the  minimum  burst-depth  conditions  must  be 
incorporated  to  finally  arrive  at  the  desired  detonation  time. 


i 

( I 
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In  determining  t,  several  possible  situations  arise.  First,  the  PCA  may  be  reaehed  at  or 
above  the  surface.  In  this  case  r s 0.  Once  the  missile  is  beneath  the  surface,  it  may  reach  the 
PCA  before  reaching  its  maximum  penetration.  In  this  ease  r = /, . On  the  other  hand,  the  missile 
might  stop  at  its  maximum  depth  before  achieving  its  PCA.  In  this  case  t = r,.  To  determine 
which  of  these  two  realistic  times.  /,  or  is  the  time  to  reach  the  closest  achievable  point,  tvvo 
critical  miss  distances,  f?,  and  are  necessary. 


The  parameters  a,  9.  t/.  A,  and  V determine  a family  of  parallel  trajectories.  The  miss 
distance  determines  a unique  member  of  a family.  A trajectory  with  miss  distance  /?,  will  be  at 
the  PCA  upon  impact.  For  the  PCA  to  be  on  the  surface,  ty  = 0.  That  is. 

V I 

- — - - — VV*  ± 2A  cos  a R,  + 2A</sin  a . (49) 

A A 

Squaring  and  simplifying, 

0 = ±2AAf,  cos  a cos  + 2A</sin  a . (50) 

Thus, 


R,  = 


d sin  a 
cos  a cos  0 ' 


(51) 


On  the  other  hand,  a trajectory  with  miss  distance  Rt  will  reach  the  PCA  at  the  maximum 
penetration  depth.  Hence  R^  will  depend  upon  the  maximum  depth  as  well  as  the  initial 
parameters.  For  the  PCA  to  be  at  the  maximum  depth,  t.  = r,.  Therefore,  the  radical  in  equation 
(46)  must  be  zero. 

V*  ± lAR^  cos  a cos  9 + 2Ad  sin  a = 0 (52) 

Solving  for  R^, 


Rx 


V*  + 2Ad  sin  al 
2A  cos  a cos  sJ 


(53) 


A particular  maxiuium  penetration  depth,  M,  forces  /?,  = 0.  From  figure  8, 

M J 

where  J is  as  shown  in  figure  8.  So  A/  = </  sin  a*. 

! V* 

The  constraint  Rf  = Q implies  sin  a = - from  equation  (52).  So 


(54) 


M = 


V* 

4A*</' 


(55) 


However,  referring  to  figure  9,  for  an  early  impact,  any  maximum  penetration,  is 

equivalent  to  ^ s 0 and  fL>  M\s  equivalent  to  fZj  > 0.  Referring  to  figure  10,  for  a late  impact, 
s Af  is  equivalent  to  ^ > 0,  and  /tt  a Af  is  equivalent  to  s 0.  Thus,  M is  not  necessary  for 
the  analysis.  Only  the  sign  of  R^  is  pertinent  to  the  analysis. 


Then,  referring  to  figure  9,  for  an  early  impact,  the  processor  calculates  the  value  of  R^ 
using  equation  (53)  with  the  upper  signs.  If  s 0,  then  all  early  impacts  will  reach  maximum 
penetration  before  the  PCA.  Thus,  t = . If  the  miss  distance  parameter  R satisfies  Q < R <■  R^, 


*b 
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the  missile  will  reach  the  PCA  before  stopping.  Then  t = r, , If  /?  a i?, , the  PCA  is  below  the 
maximum  penetration  depth.  Hence,  r = tj. 

Referring  to  figure  10  for  a late  impact,  the  data  processor  computes  the  values  of  /?,  and 
using  equations  (51)  and  (53)  with  the  lower  signs.  Three  cases  arise  when  > 0.  If  0 s ^ 

T = li;  the  missile  reaches  its  maximum  depth  before  the  PCA.  If  ^ » the  missile 

reaches  the  PCA  above  the  maximum  depth,  so  t|  = t.  If  /?  a Af,  the  PCA  is  reached  before 
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Figure  9.  Critical  miss  distances  for  early  impacts  with  various  maximum  penetration  depths. 
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Figure  10.  Critical  miss  distances  for  late  impacts  with  various  maximum  penetration  depths. 


impact  and  t s 0.  Two  cases  occur  for  f?,  ^ 0.  If  /?>/?, , the  PCA  is  again  aboveground.  But  if 
0 </?</?,,  the  PCA  is  reached  below  the  surface  but  above  the  maximum  depth.  In  this  case 

T= 

To  illt...d>e  this  behavior,  dependent  on  two  critical  miss  distances,  the  HP98I0  was 
programmed  to  plot  p*(t).  In  figure  1 1,  p*(/)  is  plotted  for  various  miss  distances  (noted  above  the 
curves)  using  the  following  values:  d - 30.5  m,  o = ^ = ir/3,  A = -3.6  x 10^  m/s*  and  V = 610 
m/s.  Assuming  an  early  impact,  the  critical  miss  distance  Rj  = 8.95  m. 
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Figure  II.  The  effect  of  R on  target  range  as  function  of  time.  : 

I 


There  is  a minimum  desired  burst  depth  constraint  to  be  considered  also.  The  time,  T,  to 
reach  this  depth,  m,  is  the  solution  to 

z(7)  = m = - -j  sin  a 7*  - Fsin  aT  + d.  (56) 

Solving  for  T gives 


r = 

The  time,  T,  must  be  the  smaller 


V Vv*  sin’  a + 2A(d  — m)  sin  a 
A ~ (-M  sin  a)  ’ 

positive  value  since  the  larger  corresponds  to  a negative  missile 


velocity.  Assuming  that  d > m,  that  is,  that  the  target  is  below  the  minimum  depth,  then 

(2A  sin  a)(d  ~ m)  <0  (58) 

since  A < 0.  Therefore,  V sin  a is  greater  than  the  radical.  So 

V sin  a + v/V?  sin*  a + 2A(d  - m)  sin  a 

r — . (59) 

A sin  a 


The  minimum  burst-depth  criterion  forces  the  logic  (see  fig.  12)  to  state  that  detonation  occurs  at 
T when  T 2 T.  Otherwise  detonation  occurs  at  T. 


INPUT 
y .K.  a.m.  A. 

V.d.  R 


compute 

(/.  « 


COMPUTE  YES  ^ 
R,.  Rj.  T, 


NO 

driMPUtE 

..  R2.  T.  tl 

• max 

R >R, Rj 


RjSo? 


H.5R2 


Figure  12.  Logic  diagram  for  linear  underground  path  to  point  target. 
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3.2  Case  11:  Linear  L'nden(ruund  Path  to  Isovulntrable  Sphere 

In  this  case,  we  consider  a linear  underground  trajectory  where  a detonation  is  desired  at  the 
PCA  to  an  isovulnerable  sphere.  The  target  is  assumed  to  be  inside  an  isovulnerable  sphere  of 
known  radius  r.  The  depth  to  the  center  of  the  sphere  is  now  d.  The  values  of  A,  V,  R,  a,  ff,  y,  ir, 
and  d are  known.  The  procedure  is  essentially  the  same  as  in  Case  I with  the  new  objective  being 
detonation  as  close  as  possible  to  the  isovulnerable  sphere.  The  distance  function  p(f)  and  the 
position  functions  jf(/),  y(/),  and  :(/)  now  are  with  respect  to  the  center  of  the  sphere.  In  addition 
to  the  logic  of  Case  1,  the  fuze  can  detonate  at  any  time  t for  which  p(r)  ^ r.  The  times  T,  and  Tj 
are  the  times  when  pit)  = r,  that  is,  when  the  missile  pierces  the  isovulnerable  sphere.  These 
values  are  the  two  real  minimal  positive  roots  of  p\i)  - /^  = 0.  But  the  solution  of  the  fourth- 
degree  equation  in  time  is  quite  involved.  The  procedure  is  outlined  in  appendix  B and  can  be 
programmed  with  impact  parameters,  if  so  desired.  The  result  is  a time  interval  [T, . Tj]  in  which 
to  detonate.  A second  method  is  for  the  fuze  data  processor  to  continually  compute  pit)  and 
detonate  shortly  after  the  time  T^ , at  which  p(T,)  = r. 

However,  for  both  Cases  I and  II,  a third  method  for  finding  the  PCA  is  to  continuously 
compute  and  compare  successive  values  of  pit).  When  the  pit)  values  stop  decreasing  and  start 
increasing,  the  missile  is  at  a PCA.  If  the  missile  is  below  the  minimum  burst  depth,  m. 
detonation  can  occur  at  this  time.  But  there  are  two  problems:  (I)  how  many  values  are  enough 
to  ascertain  that  pit)  has  begun  to  increase,  ruling  out  random  error  fluctuations,  and  (2)  how  to 
proceed  if  pit)  is  initially  increasing.  The  former  problem  requires  experimental  data  before  an 
answer  can  be  made.  The  latter  problem  could  occur  if  a is  so  small  that  the  missile  reaches  the 
PCA  before  the  data  processor  can  compute  T and  t This  difficulty  could  be  avoided  by 
providing  an  additional  condition  to  use  the  third  method  if  pit)  is  initially  decreasing.  If  not, 
detonation  should  occur  at  T. 

The  following  is  a summary  of  the  logic  for  a linear  underground  path  to  a point  target. 

1.  Input  initial  parameters  m,  y,  a.  A,  V,  a,  and  R.  If  R will  not  be  available,  input  a or  A and 

calculate  R. 

2.  Compute 

(a)  O',  e, 

ib)  form  of  impact, 

(c)  appropriate  critical  mis"  distances,  and 

(d)  T,  ti , tj. 

3.  Detonate 

(a)  at  r = T if  T > T,  or 

(b)  at  f = T if  T s r. 
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4. 


\ 

First  alternative  tn  Pt'A  fii/iiig 
(a)  Pi\)eeed  as  iii  steps  I thnuigh  2(h). 

(h)  Monitor  p(.). 

(e)  When  p(/)  reaehes  a local  tninirniini,  detonate,  if  z{t)  in. 

(d)  If  at  this  niininium  ;(/)  in.  delay  detonation  until  :(t)  ^ in. 

5.  Second  alternative,  for  lethal  sphere  ecinsiderations  f 

(a)  Monitor  p(/)  until  pit)  '-Z  r.  ' 

(b)  Detonate  when  :(i)  s i:i.  ■' 

3.3  Case  III:  Nanlinear  Underj>round  Path 

In  this  case,  the  initially  linear  trajectory  after  entry  becomes  a curve  because  of  missile 
accelerations.  As  in  the  linear  underground  trajectory  cases,  it  is  assumed  that  the  values  of  a.  y, 

(T.  fl.  R.  V,  in.  and  dare  known  before  entry.  If  the  target  is  within  an  isovulnerable  sphere,  the 
depth,  d,  becomes  ihe  depth  to  the  center  of  the  sphere.  By  hypothesis,  the  missile  has 
accelerometers  to  measure  the  ciimponents  of  the  acceleration  A with  respect  to  the  target 
oriented  coordinate  system  For  this  discu;  s on  it  will  also  be  assumed  that  the  missile  remains 
underground  after  impact.  Figure  13  illustrates  the  geometry  for  this  case.  The  parametric 
equations  of  position  wi'h  respect  tt)  time  are 


.v(t)  = f f A t.v)  d.v  dn  + Vj./  + .r„, 

•'«  •'ll 

(60) 

.v(t)  = f f ''fj/l-v)  d.v  d«  ^ V^t  + , and 

•'o  -'ll 

(61) 

;(/)  = f f ds  du  + + z,, , 

•'ll  *'0 

(62) 

/ 

r 

z 


Figure  13.  Nonlinear  underground  trajectory  to  point  target  at  center  of  isovulnerable  sphere. 
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where  (j:o*  .'o'  O*  and  the  velocity  components  are 

V,  = V cos  a cos  0, 

V,  = V sin  0 cos  a,  and 
V,  = V sin  a. 

Using  z(/)  for  depth,  the  missile  will  reach  minimum  burst  depth,  m,  when 


m = 2(71 


ds  du  - (V  sin  a)  T + d. 


(63) 


The  data  processor  must  compute  this  function  continuously  to  determine  when  the  missile  is 
below  the  minimum  depth. 


Next,  the  distance  between  a point  on  the  trajectory  and  the  target  is  given  by 


pit)  = V'x*(/)  + y*(/)  + T}{t). 


(64) 


After  substitution  this  becomes 


pit)  - J A,(5)  ds  du  + V cos  a cos  0t  + 
Ayis)  ds  du  + Vcos  a sin  6lj 
(f  f ^ti^)  ttu  - Vsinaf  + </)  j . 


(65) 


A first  approach  to  the  problem  of  minimizing  pit)  is  to  assume  that  the  path  has  only  one 
minimum.  Then  the  data  processor  should  calculate  pit)  values  continually,  as  in  the  linear  case, 
to  determine  when  pit)  is  a minimum.  If  the  missile  is  below  the  minimum  depth  at  this  point, 
then  the  warhead  should  be  detonated.  Otherwise,  detonation  should  occur  as  soon  as  z(f)  = m. 

However,  the  trajectory  may  have  several  local  minima.  If  ti  e above  procedure  is  used,' the 
resultant  minimum  may  not  be  the  closest  point  achievable.  One  alternative  would  be  to  detonate 
at  the  first  such  minimum  once  below  the  minimum  burst  depth.  Another  alternative,  which  is 
useful  if  there  is  a isovulnerable  sphere  of  known  radius,  r,  would  be  to  detonate  when  pit)  ^ r. 
A third  method,  a combination  of  the  above,  would  be  to  detonate  when  pit)  s r or  when  the 
path  reachv^s  a local  minimum,  whichever  comes  first  after  passing  the  minimum  depth. 


4.  LINEAR  TARGET  WITH  CYLINDRICAL  ISOVULNERABLE  VOLUME 

As  in  the  previous  analysis,  it  is  assumed  that  the  values  of  a,  0,  R,  y,  o-,  V,  and  m are 
known.  If  the  targe'  is  a line,  the  analysis  considers  a cylindrical  isovulnerable  contour  of  length 
/ and  radius  r,  with  hemispheres  of  radius  r on  the  ends  of  the  cylinder.  Such  a configuration  is 
shown  in  figure  14.  The  total  axial  length  is  / + 2r.  The  depth  to  the  top  of  this  volume  is  d.  The 
origin  of  the  coordinate  system  is  at  the  bottom  of  the  volume  so  that  the  aim  point  on  the 
surface  has  coordinates  (0,  0,  D),  where  D = 2r  + I + d.  The  objective  is  to  detonate  at  the  PCA 
to  this  volume.  Assuming  the  missile  is  below  the  minimum  depth,  there  are  three  possible  points 
on  the  trajectory  for  detonation,  depending  on  the  depth  of  the  missile.  TTiese  possible 
detonation  points  are 
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Figure  14.  Linear  underground  trajectory  to  isovulnerabie  volume  about  linear  target. 

(1)  thePCAto(  OJ+r), 

(2)  the  PCA  to  the  z-axis  if  the  z coordinate  of  the  missile  position  satisfies  r < z < I + r, 
and 

(3)  the  PCA  to  (0,  0,  r)  if  the  z coordinate  of  position  is  less  than  r. 

For  simplicity,  it  is  assumed  that  detonation  is  to  occur  at  the  first  such  PCA  reached  below  the 
minimum  depth. 

The  parameterizations  of  the  position  of  the  missile  are  similar  to  those  for  the  point  target. 
First,  assume  the  underground  path  is  linear.  Then, 

x(t)  is  given  by  equation  (42), 

y(t)  is  given  by  equation  (41),  and  (66) 

—A 

z(t)  = sin  - V sin  at  + D. 

Again,  the  time  to  reach  the  .ninimum  burst  depth,  T,  is  obtained  by  solving  2(7j  = m.  Thus 

— V \/ V*  sin’  a + 2A{D  — m)  sin  a 

T = — + . (67) 

A A sin  a 

To  find  the  time  to  reach  each  of  the  three  possible  PCA’s,  three  distance  functions  are  required. 
P(/)  is  defined  to  be  the  distance  from  any  point  (x,  y,  z)  to  (0,  0,  / -t  r).  M(i)  is  defined  as  the 
distance  from  any  point  to  the  portion  of  the  z-axis  between  (0,  0,  / + r)  and  (0,  0,  r).  B{t)  is 
defined  as  the  distance  from  any  point  to  the  point  (0,  0,  r). 
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The  time  to  reach  the  point  on  the  trajectory  closest  to  (0.  0,  r + /).  r, , is  a solution  to 

= (68) 

where 

/“{/)  = xHl)  + yHt)  + [2(r)  - (r  + /)?.  (68) 

After  substituting  for  x{t),  y(t),  z(/),  and  simplifying, 

/«(/)  = + AVf'  + (V**AK  cos  « cos  a -A(<^  + r)  sin  a)t^  (69) 

4 

+ 2V’(+/?cos«  cos  9 - (il  + r)  sin  a)/  + /?’  + (t/  + r)*. 

After  setting  the  derivative  equal  to  zero,  equation  (69)  is  in  the  same  form  as  equation  (44).  So 
its  solution  is  given  by 

V , ; 

T,  = — -T  VV*  + 2A  [(d  + r)  sin  a ± R cos  9 cos  a)  (70) 

A 


or 

V 

r,  = --. 

The  time  to  reach  the  point  on  the  path  closest  to  the  z axis,  r,,  is  a solution  to 

where  Af‘(/)  = jr*(/)  + y*(/).  That  is. 

At 

M*(t)  = — cos*a/^  + AV(cos*a)r’  + (V*cos*a  RA  cosa  cos  9)t*  (71) 

4 


' 2 V/?(cos  a cos  8)t  + R*. 


_After  setting  the  derivative  equal  to  zero  and  simplifying, 

n _2V7?cOS0 

0 = r’  + 3-r«  + 2 — + - t -jl • 

A LA*  A cos  aJ  A*  cos  a 

The  same  solution  technique  as  for  equation  (44)  gives  the  solutions; 


T,  ■=  - — , or  T, 


^ ^ * /■/»  ^ 1 /f  cos  9 
A A » cosa 


(72) 


(73) 


The  time  to  reach  the  point  on  the  trajectory  closest  to  the  point  (0,  0,  r),  r,,  is  a solution 
d 

to  -T  B^{t)  - 0,  where 
at 

EP(l)  = x*(/)  + y*(r)  »■  (z(r)  - r)*. 

In  teims  of  initial  parameters  B*  becomes 


P(t)  = ^r^  + AV/*  + (V*5F  RA  cos  a cos  9 - A[D  - r)  sin  a)l* 

4 


(74) 


+ 2V(-D  sin  a cos  a cos  9)t  + Cf  + R*. 
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This  expression  for  B*(/)  is  of  the  same  form  as  equation  (44)  where  d is  (D  — r)  here.  Thc 
solutions  are  given  by 


V'  I 

r,  = - r + — W*  + 2A(D  - /•)  sin  a ± R eos  a cos  S) 
A A 


or 


(75) 


The  objective  is  to  detonate  at  the  point  of  closest  approach  that  the  missile  can  achieve 
below  the  surface.  The  times  t,  . Tj,  and  are  the  times  to  the  PCA's  along  an  infinitely  long 
trajectory.  They  w !l  always  occur  in  this  order.  However,  the  missile  may  reacti  these  points 
before  'mpact.  Further,  the  missile  rray  not  reach  any  of  these  points  at  all  if  they  are  below  the 
maximum  attainable  depth.  Again  t - -VfA  is  the  time  at  which  the  missile  reaches  the 
maximum  depth.  As  in  the  poin.  target  tase,  there  are  certain  critical  miss  distance  values 
necessary  to  determine  which  of  these  times  is  optimal  for  detonation. 


Aj  and  /?4  are  the  miss  distances  for  which  the  PCA's  to  (0.0,/+r)  and  to  (O.O.r)  respectively, 
occur  at  the  maximum  penetration  depth,  f?,  and  R^  are  the  miss  distances  for  which  the  PCA’s 
to  (0,0,/+r)  and  (O.O.r),  respectively,  occur  upon  impact.  These  four  critical  miss  distances  are 
computed  i'’  the  same  way  as  their  counterparts  of  the  point  target  case.  Thus, 


R, 

R. 

R, 

Ri 


- r)  sin  tirl 

2a  cos  a cos  0 J ’ 

1 

(D  - r)sina1 
cos  0 cos  a J ’ 

+ 2A(d  + r)  sin  nj 
2.4  cos  0 cos  a J ' 

(d  + r)  sin 
cos  0 cos  a J ' 


and 


(76) 

(77) 

(78) 

(79) 


if  the  impact  is  early,  /?,  and  R^  will  both  be  negative.  These  critical  miss  distances  also 
satisfy 


Rf  Rf  (otherwise  d>  D which  contradicts  the  hypothesis). 


fi,  & /?,  (otherwise  V*  < 0),  and 


Ra  > Ri  (otherwise  V*  < 0). 

Figures  I.S  16,  and  17  diagram  the  various  relationships  for  an  eariy  impact.  The  procedure  to 
determine  ti.e  detonation  time  for  an  early  impact  is  summarized  as  follows. 


1.  If  /i,  s 0.  then  R>  Ri  for  any  early  impact  and  the  missile  stops  before  reaching  a trajectory 

PCA.  Therefore,  detonation  occurs  at  <*. 

2.  If  /^  > 0,  > 0.  and  R satisfies 

(a)  0 < /?  < /?4,  the  missile  reaches  all  trajectory  PCA's  in  order.  Detonation  occurs  at  the 
first  T a T.  If  there  is  no  t 2:  T,  detonation  occurs  at  T. 

(b)  R4  < R < Rf,  the  missile  reaches  the  ‘rajectory  PCA  to  the  top  of  the  target  first  and 
then  to  the  middle  of  the  target.  Detonation  occurs  at  whichever  t a T.  Otherwise 
detonation  occurs  at  . 
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(c)  H.  s R.  the  missile  passes  through  all  the  PCA's  of  the  path  on  or  before  impact. 
Detonation  occurs  at  T. 

3.  If  R,  > 0.  R,  s 0,  and  R satisfies 

(a)  0 < /?  < the  missile  reaches  the  PCA  with  respect  to  the  top  first  and  to  the  : axis 
last.  Detonation  occurs  at  whichever  t ^ T.  Otherwise  detonation  occurs  at 
(h)  R.  R.  the  missile  reaches  its  maximum  penetration  depth  before  any  trajectory  PCA. 
Detonation  occurs  at  /j. 


I igure  If).  Harly  impact  PCA  relationships  when  Rt  > Oand  /?<  > 0. 


V 


j 


I ' • 
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Figure  17.  Early  impact  PCA  relationships  when  /Sj  > Oand  R,  = 0. 


On  the  other  hand,  if  the  missile  is  late,  /?|  and  R^  are  both  positive.  Further,  these  miss 
distances  satisfy 

R,  2 Rj  (otherwise  J > D,  contradicting  the  hypothesis) 

R,  > (otherwise  V*  < 0),  and 
/?3  2:  Rt  (otherw  ise  V*  < 0). 

Figures  18,  19,  and  20  diagram  the  various  relationships  for  a late  impact.  The  procedure  to 
determine  the  detonation  time  for  a late  impact  is  summarized  as  follows. 

1 . If  s 0 and  R satisfies 

(a)  0 < R < R, , the  missile  reaches  all  the  trajectory  PCA's.  Detonation  occurs  at 
whichever  t > 7 occurs  first.  Otherwise  detonation  occurs  at  7. 

(b)  R,  s /?</?, , the  missile  reaches  the  PCA  to  the  top  before  impact.  After  impact  it 
reaches  the  PCA  to  the  z axis  first  and  to  the  bottom  last.  Detonation  occurs  at  t*  or  Tj, 
whichever  is  greater  than  7.  Otherwise  detonation  occurs  at  7. 

(c)  Rt  < R,  all  the  PCA's  of  the  path  occur  before  impact.  Detonation  occurs  at  T, 

2.  If  > 0,  s 0.  and  R satisfies 

(a)  /?</?,,  the  missile  reaches  first  the  PCA  to  the  top,  then  to  the  z axis.  Detonation 
occurs  at  t,  , or  t,  , whichever  is  greater  than  7.  If  neither  is  greater  than  7,  detonation 
occurs  at  I2 . 

(b)  /?,  ^ R<  R,,  the  PCA  to  the  top  occurs  before  impact.  After  entry,  only  the  PCA  to  the 
z axis  occurs  before  the  missile  stops.  Detonation  occurs  at  Tj  if  T2  2 7.  If  not, 
detonation  occurs  at  . 

(c)  Rt^  R ^ R^,  the  missile  after  entry  reaches  the  PCA  to  the  z axis,  followed  by  the  PCA 


\ 
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Figure  20.  Late  impact  PCA  relationships  when  > 0 and  /Zj  > 0. 


to  the  bottom.  Detonation  occurs  at  r,  or  Tj,  whichever  is  greater  than  T.  Otherwise, 
detonation  occurs  at  . 

(d)  R3  < R,  all  the  PCA's  of  the  path  occur  before  impact.  Hence  detonation  occurs  at  T. 

3.  If  > 0,  /Z,  > 0,  and  R satisfies 

(a)  0 < /Z  £ the  missile  reaches  its  maximum  depth  before  any  PCA.  Detonation  occurs 
then  at  r, . 

(b)  yZj  £ yz  £ yZ,,  the  situation  is  the  same  as  2(a)  for  a late  impact.  | 

At  this  point,  this  case  becomes  identical  to  (2).  Steps  3(c),  (d),  and  (e)  correspond  to  2(b),  (c), 
and  (d).  The  complete  logic  flow  chart  for  a linear  path  to  a linear  target  is  shown  in  figure  21. 

With  all  other  assumptions  holding,  now  the  underground  path  is  taken  to  be  some  9urve 
determined  by  the  independent  interaction  of  nonuniform  component  accelerations  in  the  ^get 
coordinate  system  as  shown  in  figure  22.  Then  the  x(t),  y{t),  and  z(r)  functions  are  found  by 
integration  of  these  accelerations.  So  x(r),  y(/),  and  z(/)  are  as  defined  in  equations  (60),  (61),land 
(62),  where 

V,  = V cos  a cos  6,  and  = R, 

V,  = V cos  a sin  0,  y©  = 0, 

V;  = Vsina.  z;,  = D = / + 2r  + d. 

To  determine  th-;  depth  of  the  missile,  the  data  processor  must  evaluate  z(/).  The  minimum  depth 
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is  attained  when 


z{T)  - m - f f ils  dll  - V(sin  T + D - m = 0. 

•'ll  •»n 


(80) 


The  distance  functions  M'Ht),  /“(/),  and  BHi)  are  defined  as  before,  but  they  are  now  expressed 
with  the  new  forms  of  .r(r),  .v(/),  and  z(t).  For  example. 


mt) 


= [ f f Aj(s)  ds  dti  -t-  ^^(cos  « cos  0)t  + /?! 
L Jg  Jg  J 

^[ff  /4^{.v)  ds  du  + V(cos  « sin  0)/j  . 


(81) 


There  are  two  possible  methods  of  determining  the  point  for  detonation  but  both  entail 
monitoring  the  distance  functions  and  z(t),  as  outlined  below. 


I.  Method  A. 

(a)  If  /"(/),  and  B’(t)  are  all  increasing  initially,  detonate  at  T. 

(b)  If  r + / < ;(/)  < m and  \P(t)  is  increasing,  detonate  when  /“(t)  ^ or  /"(/)  changes  from 
a decreasing  to  an  increasing  function,  whichever  occurs  first  below  the  minimum  depth. 

(c)  If  r s zU)  i + r,  detonate  when  Af^fr)  s or  when  Af'(r)  changes  from  a decreasing  to 
an  increasing  function,  whichever  occurs  first  below  the  minimum  depth. 


Figure  22.  Nonlinear  underground  trajectory  to  isovulnerable  volume  about  linear  target. 
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(d)  If  ;(/)  < /•.  detonate  when  BHt)  ^ r*  or  when  B*{t)  changes  from  a decreasing  to  an 
increasing  function,  whichever  occurs  first  below  the  minimum  depth. 

(e)  Otherwise  detonate  when  :(t)  = m. 

2.  Method  B. 

(a)  If  ff'(f).  and  /”(/)  are  all  initially  increasing,  detonate  at  T. 

(b)  If  one  of  these  three  distance  functions  is  initially  decreasing,  detonate  when  that 
function  starts  increasing,  provided  :(/)  ^ m. 

(c)  If  more  than  one  of  these  distance  functions  is  decreasing  initially,  detonate  once  below 
the  minimum  depth  at  the  first  time  for  a function  to  change  from  decreasing  to  increasing. 

5.  CONCLUSIONS 

The  assumed  objective  of  the  fuzing  algorithm  is  to  detonate  the  warhead  at  the  closest 
achievable  point  to  the  target  when  beneath  the  minimum  burst  depth.  In  the  analysis  it  is 
assumed  that  the  missile  remains  underground  after  entry.  In  most  of  the  discussion  all  the 
necessary  parameters  are  assumed  to  be  known  precisely.  Terms  that  appear  to  be  undefined  as 
a parameter  approaches  zero  are  well  defined  if  they  are  expanded  as  a series.  The  fuze 
algorithm  results  obtained  for  the  assumption  of  only  a known  distribution  of  miss  distance  are 
briefly  discussed  in  appendix  B. 

The  general  logic  for  the  fuze  in  all  cases  is 

(a)  receive  the  impact  data  (a.  y,  cr,  f',  ^4,  R,  and  d), 

(b)  determine  which  distance  functions  are  appropriate, 

(c)  compute  and  choose  the  applicable  time  t to  reach  a PCA, 

(d)  compute  the  time,  T,  to  reach  minimum  depth, 

(e)  compare  the  values  of  T and  r,  and  then 

(f)  detonate  at  whichever  time  is  the  greater. 

The  logic  for  the  linear  path  to  a target  can  be  analytically  calculated  for  a deterministic  answer. 
For  the  nonlinear  case,  the  fuze  must  repetitively  calculate  various  distance  functions, 
comparing  successive  values  with  previous  values  to  determine  when  the  missile  reaches  a PCA 
or  the  minimum  depth,  or  when  the  missile  pierces  the  lethal  volume. 

For  a complete  model,  much  more  must  be  incorporated  into  it.  The  fuzing  logic  should 
include  a'lowances  for  (1)  the  prediction  of  a subsurface  path  from  known  soil  characteristics,  (2) 
a trajectory  that  oscillates  above  and  below  the  ground,  (3)  nonsymmetric  targets,  and  (4) 
inaccurate  impact  data.  Further,  the  fuzing  program  must  be  flexible  enough  to  be  applicable  to 
any  form  of  target  under  the  varied  conditions  of  impact,  soil,  and  data  availability. 
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APPENDIX  A 

DIAGRAMS  OF  PARAMETER  ANGLE  RELATIONSHIPS' 

In  the  computations  for  the  fuzing  logic,  the  surface  plane  angle  ^is  required.  Upon  impact, 
though,  only  available  are  the  map  cc ordinate  angles  y.  measured  countc.clockwise  from  north 
to  vector  miss  distance  R.  and  <r,  measured  counterclockwise  from  north  to  S.  the  surface 
component  of  the  trajectory.  To  formulate  an  expression  for  0 using  y and  tr.  as  intermediate 
parameter  is  required.  The  angle  O'  is  defined  to  be  the  smaller  angle  between  R and  S. 
Further.  O'  determines  whether  the  impact  is  early  or  late. 

Figure  A-t(a)  through  (x)  represents  the  surface  map  coordinate  system  at  the  point  of 
impact.  In  each  diagram  a possible  location  of  the  aim  point  is  plotted  together  with  a vector 
representing  a possible  surface  component  of  the  trajectory.  The  map  angles  yand  (rare  shown, 
as  is  *he  desired  parameter  O'.  Hence,  figure  A- 1(a)  through  (x)  represents  all  possible  relations 
between  <t,  y,  and  O'.  In  each  section  of  this  illustration.  O'  is  obtained  numerically  in  one  of  two 
ways. 

If  0 ^ |<r  - y|  s TT,  then  ^'  = |<r  — y|. 

Ifff  s |o- - y|  s 2fl-,  then  O' = 2w  - ftr  - yj. 

In  all  cases,  0 s s rr. 

Each  of  these  figures,  when  analyzed  as  to  the  type  of  impact,  early  or  late,  yields  criteria 
dependent  on  O'  for  the  type  of  impact.  The  results  are  shown  in  tables  A-l  and  A-II. 


TABLE  A-I 

Dependence  of  Impact  Form  on  O'  when  S and  R are  in  Nonadjacent  Quadrants 
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Figure  A-I.  Initial  parameter  angle  relationships  (cont'd). 
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TABLE  A-» 

Dependence  of  Impact  Form  on  O'  when  S and  R arc  in  Adjacent  Quadrants 


5,  APt  Early/Late 


0 

0 

L 

= 

2 

3 

L 

< 

1 

L 

< 

4 

L 

< 

L 

< 

3 

L 

< 

3 

0 

E 

> 

4 

L 

< 

1 

E 

> 

2 

E 

> 

1 

0 

E 

> 

3 

L 

= 

1 

L 

= 

4 

L 

< 

2 

L 

< 

3 

L 

< 

4 

0 

E 

> 

4 

L 

< 

I 

E 

> 

2 

E 

> 

2 

0 

E 

> 

3 

E 

> 

1 

E 

5 • 

4 

L 

2 

L 

ss 

Figure  A-2,  a composite  generalization  of  figure  A-l(mJ  through  (t),  shows  the  possible 
relationships  between  S and  R when  in  acjacent  quadrants.  Ri  is  the  vector  miss  distance  from 
. the  /th  aim  point,  AP, , to  the  impact  point.  Sf  is  the  surface  component  of  the  ith  trajectory.  Any 

relationship  shown  in  figure  A- I(m)  through  (t)  is  represented  by  a pair  (5( , Pj)  for  some  (/,»  pair 
^ where  /,  y = 0,  1,2,  3,  4. 
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APPENDIX  B 

CALCULAIION  OF  EXPECTATION  FOR  MISS  DISTANCE 


If  {he  miss  distance  R is  not  known  precisely,  the  values  for  R arc  assumed  to  be  Rayleigh 
di'^tributed  with  standard  deviation  rrand  distribution  function  F(R)  given  by 

F(R)  = (B-l) 

tr 

The  analysis  seeks  the  time  to  reach  the  point  of  closest  approach.  Since  /,  depends  on  R.  if  R 
is  not  known  precisely,  the  expected  value  of  R,  E{R).  is  used  instead  to  calculate  . Using  the 
variables  V’for  the  impact  velocity,  dfor  the  depth  to  the  target.  A for  the  acceleration.  « for  the 
impact  angle,  and  rtfor  the  surface  plane  angle. 

-V  1 f* 

£(/?)  = ——  + — I [V'*  + 2/t(±/?cos  fleosa  + </sin  a)’'"  £(R)  (//?.  (B-2) 

A A J„ 

where  (for  ± and  + signs'  the  upper  sign  corresponds  to  an  early  impact  and  the  lower  to  a late 
impact.  Let  .r  = R and  consider  the  integral 

j VV-’  X IxA  cos  o cos  D + 2AJsm  a j </.f.  (8-3) 


The  data  processor  can  estimate  this  integral  by  numerical  methods.  An  alternative  is  to  estimate 
(B-3)  by 

f dx.  (B-4) 

•'O 

Since  the  order  of  magnitude  of  V*  is  I0^  while  the  other  terms  in  the  radicand  contribute  10*. 
the  radical  is  approximated  by  V.  But 

I ^ </jr  = V.  (B-5) 

~V  V 

So  E(R)  is  approximately  — r-  + — = 0.  The  other  times  computed  for  the  various  cases  all  have 

A A 

this  ^orm.  Therefore  the  expectation  of  each  time  will  be  zero,  suggesting  that  this  is  too  simple 
an  approximation  to  equation  (B-3).  Numerical  approximation,  say  by  Simpson's  Rule,  appears 
to  be  a better  method. 
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APPKNDIX  C 

SOLLH  ION  OF  A FOUR  f II-DFGRFF  KQUATION 


F or  a liiu'ar  uinJcii!roiiiKl  path  to  a Ictli;,!  sphere  of  radius  r.  the  object  is  to  ascertain 
x'.hen  the  missile  is  vsithin  the  lethal  sphere.  Hence  the  limes.  1.  for  which  frit)  < arc 
necessary  where  pit)  is  the  distance  to  the  center  of  the  sphere  as  a function  of  time.  To  find 
such  values  of  i.  the  equation  frit)  - r = 0 must  be  solved  for  t.  Using  the  variables  V for 
imrael  velocity,  J for  depth  to  the  target.  A for  acceleration,  a for  impact  angle,  and  0 for  the 
suitace  plane  angle,  frit)  ~ r - 0 bcreomes  as  follows. 


0 = t'  + 
+ X 


4V1  ru- 

— + 4 — 

A i L /f  - 

V 

— iR  cos  H cos  a + (/  sin  o) 


R cos  0 cos  o d sin  i 


4 

t + — iR- 
A" 


+ t/-'  - r). 


(C-l) 


where  (for  e.  ;ind  ■+•  signs)  the  upper  sign  corresponds  to  an  early  impact,  and  the  lower  to  a late 
impact.  This  equation  then  is  of  the  form 


0 = /^  + at'  + ht-  + ct  + / 


(C-2) 


where  a.  h.  c,  ;ind  / are  arbiintry  constants.  Tv)  solve  this  equation,  first  the  following  cubic 
equtition  is  solved 

y ' - /;v'^  + iac  - 4/jy  - d’f  + 4hf  - < - = 0 where  a.  h.  c./are  as  in  equation  (C-2)  (C-3) 

Tv)  solve  (C—!).  let 

p = -h 


<1  = ac  - 4f 
r = —d-f+  4hf  - 


Define 


(C-4) 


r = - Oq  - ir) 


/ 

R = — (2/)  ’ - 9pq  -t- Hr) 

’ = K 

+ L — p/?'  where 

R 

+ 

2 

B 

4 27 


Using  this  solution  for  y.  determine  R by  evaluating 


If  R + 0.  then 


and 


R 


h + V 


R-  -2h  + 


'4ah  - Xc  - ir' 

I Jr 


ha-  j 


4ah  — Xc  — 


4R 


(C-5) 


(C-6) 


(C-7) 


(C-8a) 


(C-8b) 


Trui:3  PAG3  BlJ/'WfOT  HU 
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APPENDIX  C 

If  /?  = 0.  then 


and 


F=yl^-2h  + 2\/y^^^- 

/3«2  , 

E = yj—  -2h-2 


{C-9a) 


(C-9b) 


Then  the  solutions  of  f^U)  - = 0 are  given  by 

4 2 “ 2 


(C-10) 


The  two  real  minimal  positive  values  resulting,  call  them  T, . T^,  represent  the  times  at  which  the 
missile  trajectory  enters  and  exits  the  isovulnerable  volume.  Thus,  there  is  a time  inter  al  fT, . 
Tj]  in  which  to  detonate  effectively. 
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